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Abstract 



In a medium composed of scalar particles with non-zero mass, the range of Van-der- 
Waals-type scalar mediated interactions among nucleons becomes infinite when the 
medium makes a transition to a Bose-Einstein condensed phase. We explore this 
phenomenon in an astrophysical context. Namely, we study the effect of a scalar 
dark matter background on the equilibrium of degenerate stars. In particular we 
focus on white dwarfs and the changes induced in their masses and in their radii. 



Dark matter is one of the outstanding problems in modern cosmology 1 . Evidence for it 
comes from a multiplicity of sources and distance scales [21 E] . The nature and composition 
of dark matter might be diverse. It might well be that each different distance scale has its 
predominant type of dark matter. Several candidates for dark matter have been postulated 
jlj. The list includes the lightest supersymmetric particle (e.g., the neutralino) and the axion. 
They might be heavy, like the neutralino or light, like the axion. They might constitute the 
bulk of dark matter, like WIMPS (weakly interacting massive particles) in CDM (cold dark 
matter) scenarios or be just a small fraction of it, like neutrinos (they have mass and thus they 
constitute a fraction of the mass in the universe). Because scalar particles, such as the axion 
or the Higgs boson, are fundamental ingredients of the standard model (SM) of elementary 
particle physics and completions thereof, scalars might constitute also a natural ingredient of 
the dark stuff in the universe. 

The characteristic feature of dark matter in general is precisely the fact that it only 
interacts weakly with ordinary matter and therefore manifests itself mainly gravitationally. 
So far, evidence of dark matter is indirect. Nevertheless, direct searches that rely on the 
(feeble) interaction of dark particles such as WIMPS with ordinary matter are under way and 
there is hope that sufficiently sensitive experiments will eventually lead to positive results 
0. But there might be other effects of the feeble interactions of dark matter with ordinary 
matter that show up in other settings. In this paper we explore the consequences of scalar 
dark matter in the astrophysical environment of compact stellar objects. 

Suppose very generally that in a grander unified scheme beyond the SM, new scalar par- 
ticles carrying a new conserved quantum number Q do exist. If, at some high energy scale, 
those putative scalars share fundamental interactions with ordinary quarks and with new fields 
associated to the high energy scale, then, at low energies, the 2-body t-channel-exchange of 
these scalars will generate feeble residual spin independent forces (i.e. Van der Waals type 
dispersion forces) [HI Ej among nucleons that add coherently over unpolarized bulk matter and 
extend over distances on the order of the Compton wavelength of the mediating scalars. Fur- 
thermore, if matter is embedded in a Bose-Einstein condensate composed of these scalars, 
then the dispersion force becomes infinitely ranged, i.e. it extends over far larger distances 
than the Compton wavelength. This phenomenon was first discovered in [Hj and exploited in 
in a cosmological setting. 

Our starting point is the assumption that galactic dark matter is composed at least in 
part by (light) scalar particles. We focus our attention in a region where ordinary matter has 
accreted and evolved into a white dwarf and the condensed scalars have been gravitationally 
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trapped by this stellar matter. So, we start with an ideal boson gas in a gravitational field. 
The gas is at a temperature T ^> m where m is the mass of the scalar. The gas is relativistic. 
The statistical mechanics of a free ideal boson gas has been given in [H]. Here, we shall 
describe the statistical mechanics of Bose-Einstein condensation of a relativistic boson gas in 
an external field U(x, y, z). For condensation to occur, a non zero chemical potential \x for the 
gas is mandatory. We assumed before that our scalars carry a conserved quantum number Q. 
Particles and antiparticles, i.e. both Q = +1 and Q — — 1 scalars, constitute the gas. The 
net total charge Q is then 1 
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exp (5{E - n) - 1 exp (3{E + //) - 1 
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with E = \Jp 2 + m 2 + U(x, y, z) where U(x, y, z) ~> and vanishes at the origin, and j3 = T 1 . 
In the continuum limit Q reads 
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dV / d 6 p 



(2) 



exp p(E - (i) - 1 exp @(E + fj) - 1 _ 

In a central potential U(r), the case we shall entertain to keep things simple, this can be cast 
in the following form, 
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Q = hj £ (E)dE 
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_expP(E — fi) - 1 exp (3(E + - 1 

with the density of states e(E) = J U r 2 dr(E — U(r))y/(E — U(r)) 2 — m 2 and U(E) is such 
that U(U) = E. 



Above a critical temperature T c defined implicitly by the equation 

Q = Q(fi = m,T = T c ) 



(4) 



the charge is correctly given by equation (2). But for T below T c , equation (2) gives only the 
charge distributed in excited states. For T < T c , a macroscopic fraction of the total charge 
sits in the ground state and is not included in the above integral. The boson gas is condensed, 
and the particles in the zero momentum mode constitute the condensate. For an harmonic 
oscillator potential (i.e. U(r) = ar 2 ), as would be the case of a scalar particle inside a uniform 
matter distribution, the condensation temperature can be easily calculated to be 
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up to corrections 0(m/T). 



1 For defmiteness we take Q > 0, i.e. particles outnumber antiparticles. 
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The condensed charge in this case is given by 



Q = Q l- 




T \ 2 



(6) 



Although Bose-Einstein condensation is condensation in momentum-space, in the presence of 
a gravitational field there will be condensation in ordinary space as well, and two separate 
phases will emerge just as in ordinary gas-liquid condensation. This fact causes the charge 
Qo to accumulate at r = 0. The thermal particles, i.e. those distributed in excited states, 
bearing a kinetic energy on average large compared to their rest mass will eventually escape 
from the star. Only the condensate will be trapped at the center. Because all the charge Qo 
sits at r = 0, the density of charge grows to infinity. This is however not realistic since the 
uncertainty principle prevents both momentum and position of the condensed particles to be 
simultaneously sharp. In the gravitational field of the star, the condensate will occupy a finite 
region centered around r = whose size depends on the depth and width of the potential 
well, and on the mass m of the scalar particles. The region will be a macroscopic fraction 
of the star's volume for certain ranges of these parameters (see below). Thus, matter in this 
region will be subject to the pull of macroscopic 2-scalar exchange forces on top of the pull of 
gravity. Let us give a succinct briefing on dispersion forces in a thermal bath of the mediating 
scalars. The derivation closely follows [Hj. Consider two nucleons a distance r apart and at 
rest relative to the scalar heat reservoir. We describe the effective interaction of nucleons with 
scalars with the lagrangean 



where is the nucleon mass (chosen to set the scale of the coefficient of the dimension— 5 
effective low energy operator) and ip and are the nucleon and scalar fields, respectively. It 
is the simplest low energy realization of the underlying fundamental interactions (see [5] for a 
discussion). The force between the two nucleons arises in this approximation from the emission 
at a time by one nucleon and absorbtion by the other of two scalar quanta. We use real time 
finite temperature field theory to calculate the corresponding Feynman amplitude (see e.g. 
[TT1] for an early application of the technique). The Fourier transform of its non-relativistic 
limit is the potential we seek. 

The physical clue to the nontrivial influence of the particles in the background on the 
force due to 2-scalar t-channel exchange is that while one and two real particles cannot be 
exchanged in the t-channel, the exchange of one real plus one virtual particle in the t-channel 
is kinematically allowed. In the case at hand, the heat bath supplies the real particle. This 
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mechanism contributes a piece to the potential that reads 

1 



V ( r ) = 92 1 1 fdV r p2dp / 27r ' 2 sin2pr 



16nm 2 N rV J J yjj? + m 2 2pr 
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exp p(E - fi) - 1 expP(E + fi)-l 



Below T c a macroscopic fraction of the charge carried by particles in the reservoir piles up 
in the zero mode state (the condensate) and the integral over states in equation (8) just as 
in equations (2) and (3) no longer correctly describes the physical situation. To evaluate the 
condensate contribution to the potential in the degenerate case (T < T c ) we next establish 
its relation to Q. Call n±(p) the following distribution functions, 

1 1 
exp p(E — fx) — 1 exp p(E + fi) — 1 

Then, for T < T c , it holds that 

n+ (0) = n_(0) + — — — . (10) 

The condensate contribution to the potential is proportional to h J dVn + (0). Since U > we 
see from the previous relation, equation (10), that this quantity differs from the charge Qq in 
the condensate by less than 2(exp2/3m — 1)~ . Furthermore, (3m 1, and therefore, as long 
as the net charge Q is a macroscopically large number many orders of magnitude larger than 
T/m, the factor ~ J dVn + (0) coincides essentially with the condensate contribution to the 
charge Qq. As a consequence we can write the potential contributed by the zero momentum 
modes as (the passage from discrete to continuum momentum space, or viceversa, entails the 
substitution rule V~ l <->• (2n)~ 3 d 3 p) : 

V cond {r) = (11) 



16^777,^ m r 



where q = and the m in the denominator comes from (p 2 + m 2 ) 2 (2pr) 1 sin 2pr 
In the specific case at hand, using equations (5) and (6), we find 
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So for a degenerate scalar background the resulting potential behaves as r -1 , i.e. it is a 
potential of infinite range. On the other hand, when the temperature of the background is 
above T c , the potential is cut off by the typical Yukawa damping factor exp(— 2mr) as can be 
seen from performing the momentum integral in equation (8). 

Finally, in a vacuum, the potential takes the form jH] 

v-M = -i^V^ 1 (13) 



with Ki(2mr), the modified Bessel function whose asymptotic form gives rise to the charac- 
teristic exponential damping. We will use this result later on. 

Now that we have the form of the potential associated to the 2— scalar force among nucleons 
in a scalar heat bath and in vacuum, we turn to our problem, namely to stars whose central 
core (or even their whole volume) is permeated by a fluid made of scalar particles. To be 
definite and simple, we shall discuss equilibrium configurations called polytropes, i.e. with 
the equation of state of the form P = Kp T . More specifically, we shall analyze idealized 
white dwarfs where the Fermi pressure of a degenerate ideal gas of ultrarelativistic electrons 
supports these stars against collapse (a discussion on ordinary white dwarfs can be found in 
|12j). In the core of actual white dwarfs the nucleons are not free but bound in nuclei. In 
turn, these nuclei 2 tend to form crystal like structures, i.e. they arrange themselves on a 
regular lattice. Because of this circumstance the ideal Fermi gas equation of state is modified 
in real life by the electrostatic interactions among the positive ions and the electrons. The 
resulting corrections can be accounted for (see e.g. ^2]) and are generally small. Another 
source of departure from the ideal gas approximation is provided by inverse (3— decay which 
tends to modify the number density of the electrons and thus lessens the Fermi pressure of the 
electron gas. Here we are not concerned with this correction nor with the Coulomb correction. 
As far as the putative long range forces that we contemplate, the fact that nucleons (inside 
nuclei) cluster on lattice sites is irrelevant for our purposes since the potential associated 
to macroscopic bulk matter results from the coherent superposition of the underlying long 
distance potential created by the individual nucleons as, incidentally, it is the case for gravity. 
Therefore, to study bulk properties, the matter density is averaged over scales much larger 
than the inter-ion spacing and it is smoothed out over a macroscopically large number of 
nucleons. For neutron stars, the actual equation of state deviates from an equation of state 
for an ideal gas in ways that are more dramatic and much less under control than it is the case 
for white dwarfs. As the authors in ^2] put it, while for white dwarfs, observations of masses 
and radii are used as confirmation of astrophysical models, for neutron stars observations of 
masses and radii are instead used to test theories of nuclear physics. It is in this spirit that 
we keep things simple and do work with white dwarfs only. 

Our aim is to study the effect on stellar equilibrium of the condensation of the dark back- 
ground particles. The hydrostatic equilibrium in our case is described by the two equations: 




(14) 



2 Typically, carbon, oxygen, i.e. not exceeding atomic weight about 20. 
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for < r < L, and 

for L < r < R. 

In equation (14), 



( r l^L\ = _ 47rGp 

r 2 dr \ p dr J ' 



(15) 



K = 



9' 



(16) 



and G in both equations is Newton's constant. The distance L sets the size of the region 
where the particles of the condensate are confined. Inside a sphere of radius L, the inverse 
square law forces of gravity and 2-scalar exchange are operative. Outside this sphere only the 
pull of gravity balances the Fermi pressure. The boundary conditions for equation (14) are 
p(r = 0) = p c (where p c is the central density of the star) and p'ir = 0) = 0. The boundary 
conditions for equation (15) follow from the requirement that p and p' in equation (15) should 
match p and p' of equation (14) at r = L. The radius R of the star is implicitly defined by 
the condition p(R) = (and hence, P(R) = 0). 

_ 1 2 

(/i e is the mean molecular weight per 



In the case under scrutiny, r = | and K 



1 2 
3 — 3" 



electron) [12] ■ By changing variables to dimensionless quantities £ and 6 defined as 



p = Pc e 3 



(17) 



where the length parameter a 
in the form: 
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ttG{1+k)p£ 
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we can put both equations, (14) and (15), 



2 dX\ 2 



X 
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de 

dXi ; 



-0 ? ' 



(18) 



where X\ = £ for < £ < - and x 2 = £/vl + « for - < £ < ^. The boundary conditions in 
this new form now read: 

= (19) 



6{ Xl = 0) = 1; — 



xi=0 



for the first equation, and 



d6 

6(x 2 = L/ay/i + k) = B{x x = L/a); — 

dx 2 



de 



= -r- VI + k 

x 2 =L/a^T+K dX l 



(20) 



for the second equation. 

These equations are to be solved numerically and we shall do so below. We shall obtain the 
radius and the mass of equilibrium configurations as a function of the size of the condensate, 
i.e. as a function of the parameter = L/a. Before doing so let us recall that in the 
ordinary case, i.e. gravity alone at work, one gets, in the polytropic approximation, the 
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Chandrasekhar limit or maximum possible mass of a white dwarf. The explicit expressions 
for the Chandrasekhar mass and radius are, respectively 



Mi 



and 



where £1 satisfies = 0. 



KX* -| 

R-Ch — [ —q I Pc £1 



(21) 



(22) 



The numerical integration of the corresponding differential equation gives in this case 
& = 6.89685 and (&) | = 2.01824 12 . When plugged in equation (21), one obtains the 
well known value 

M Cft = 1.457 M . (23) 



It will be convenient for our analysis to normalize radii and masses relative to Ren and Mch, 
respectively. Trivially, 

R T„ r, 

(24) 



where X2,r verifies 6(x2,r) 



Rch 6. 

0. As to the mass, it its easy to show starting from M 



f 4:7ir 2 pdr that 



M 



M, 



2.01824" 1 
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de 




X 2,R 


dx 2 



K 



:i+S Xl ' L 



dd 




dxi 





(25) 



with xi t L = £l = ^- Sticking to ratios -j^- and make our results less dependent on the ap- 
proximations made and make them more reliable when trying to extract general consequences 
for realistic stars. 

In order to proceed we first should estimate the size of the stellar volume that can actually 
be filled by the particles of the condensate. For that purpose let us study the quantum 
mechanics of a scalar particle occupying the ground state in the gravitational potential well 
of the star. 

First of all we give the form of the potential $. For a spherical symmetric mass distribution, 

Gm(r) 



+ G / Aurpdr 



(26) 



where m(r) is the mass inside r and $(0) = $'(0) = 0. In our case, it can be explicitly cast 
in terms of the solution 9 of equations (18), 



$ = 47iGa 2 p c (l - 



(27) 
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inside the star, i.e. for r < R, while for r > R, 

$ = AitGa 2 p c + MG - ^ . (28) 

Since we shall deal only with the lowest energy mode, the wavefunction \l/(r) is purely radial 
with no orbital part. The radial Schrodinger equation can be put in an equivalent one- 
dimensional form with no centrifugal barrier, if we trade w(r) = r^(r) for ^(r): 

1 d 2 



2m dr 2 



+ m$(r) u(r) = Eu(r). (29) 



Because of the finiteness of \l/(r), u(r) should vanish at r = 0. Instead of using the actual 
$(r) (given by equations (27) and (28)), for which we have no analytical expression, we shall 
introduce in equation (29) a potential of the form <3>(r) = 7tanh 2 |. The parameters 5 and 7, 
width and depth of the well, can be fitted to give a fair approximation to the real potential. Of 
course, 5 ~ O (R) and 7 ~ O ($(00)) with $(00) being the asymptotic value of equation (28). 
The reason for doing so is that the Schrodinger equation above can be given an exact analytic 
solution. Indeed, the physically meaningful ground state wavefunction that complies with 
the boundary conditions at r = and at infinity turns out to be, after some mathematical 
manipulations, 

\U(r)=A( r -) sinh T - cosh~ 2A T - (30) 
\o/ 

-1 + y/l + 8mV 2 ] . 



with A a normalization constant and A = 7 



4 

state energy reads: 

E = (4m5 2 )' 1 -5 + 3^1 + 8m V 2 ] • (31) 

For 8m 2 7<5 2 ^> 1, equations (30) and (31) approach closely the wavefunction and the ground 
state energy of an isotropic harmonic oscillator of angular frequency u = In general, 

however, the ground state lies lower in energy and the wavefunction is broader than for the 
oscillator well. Thus, the position uncertainty of the harmonic oscillator is an underestimate 
of the uncertainty in position of the particle in the actual potential well. We shall use this 
conservative quantity in our estimates. The precise meaning we give to the parameter L is 
the following: 

with Ax 2 = (x 2 ) — (xi) 2 and r = (x±, x 2 , x 3 ). For an harmonic oscillator, Ax 2 = (2mcj) _1 . In 
our case, u = Hence, 

' ^ (33) 



n 3 1 1 
/474m2 
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L.7X 10- (1 + «)» (* \ (#- V (34) 



and, 

r, ~ 7 v m~ 4 n -u «-u _ . , , , 

\10 8 g/cm 3 / \ m J \Rch, 

if we set roughly 5 ~ R and 7 ~ 4nGa 2 p c . 

Clearly, appreciable L's —relative to stellar dimensions— will be attained only for suffi- 
ciently large k and/or sufficiently low m. However, neither can we go arbitrarily low in mass 
if we want to keep the momentum uncertainty of the scalars to be much less than their rest 
mass, nor can we increase the strength of the potential without eventually inducing unwanted 
effects in laboratory Cavendish-type experiments. We shall address these points later. Be- 
cause L depends on R, the boundary conditions (20) for the matching of the hydrodynamic 
equilibrium equations depend on the solution of those equations. Hence, if one solves equa- 
tions (18) as a function of the free parameter L and then obtains R(L) through equation 
(24), a stellar configuration should have a radius R* such that the self-consistency condition 
is satisfied: 

R* = R(L = L(RJ) (35) 

where L(R) stands for equation (34). 

We shall now integrate numerically equations (18). Our strategy is the following. We 
choose a large value for k. We do the numerics with this k fixed and letting the parameter 

, to be introduced 



de 



dx2 



and 

X 2,R 



,10 



£l vary freely. The output are the quantities X2,r, 

x 2, R ~~~ X 1,L 

in equations (24) and (25). We repeat this procedure for various values of k. To comply 
with the self-consistency requirement in equation (35) we check a posteriori for each value of 
£l = L/a whether a physically reasonable mass m can be found such that equation (34) is 
satisfied. Our results will be then presented in the parametric form a^"(^ l )}- 

It follows from the analysis that our stars tend to fall in two qualitatively distinct cat- 
egories. Those whose mass is somewhat below the Chandrasekhar mass but their radius is 
larger than the Chandrasekhar radius constitute a first group. In the second category we 
find stars with very small masses and radii. In the first case, the star adjusts itself such that 
only a small core of the star is filled by the scalar condensate whereas in the second case, 
equilibrium is attained with a large fraction of the star occupied by the condensate at the 
expense of a drastic reduction in mass and radius. The first class corresponds to low values 
of the parameter £ L (roughly up to £ L 2) and the second class corresponds to the higher 
segment of the £ L range (up to £1 = 6.89685 ; recall that £ L varies from to £1). For the 
purpose of illustration we present the two domains in two separate (R, M)-plots. In figure 
1 we show our results for 1 + k = 10 2 , 10 3 , 10 4 , 10 5 and 10 6 as a function of the parameter 
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5 10 15 20 

R/ Rch 



Figure 1: From top to bottom: equilibrium configurations for 1 + re = 10 6 , 10 5 , 10 4 , fO 3 and 
10 2 , respectively. Along the curves, £l grows from at point (1, 1) to 2 at the opposite end 
of the curves. 

< £l < 2. In figure 2 we display the results as a function of 3 < £l < £i. In this second 
plot we draw a single curve. This is because it turns out that, for large re and > 3, the 

1 3 

radius and the mass scale as (1 + re)~2 and (1 + re)~2 , respectively (this fact is reflected in the 
units chosen for the axes). Each point on the curves in figures 1 and 2 gives the radius and 
maximum possible mass of a white dwarf when the scalar condensate extends over a region 
in the stellar interior characterized by the value of the parameter £ L on that particular point 
in the curve. In figure 1, all curves merge at point (1, 1) as they should since this corresponds 
to the limit £l — > for which no macroscopic fraction of the star's volume is occupied by 
the condensate and thus the extra dispersion force is not operative. Figure 2, in turn, shows 
the other limit (i.e. £l — > £i! point (1,1) on the top left corner) for which the whole star 
is subject to the pull of the new force (gravity being negligible). A situation that is exactly 
equivalent to the Chandrasekhar case but with a blown up Newton constant. All equilibrium 
configurations on the curves of figure 1 (beyond R ~ l.bRch) correspond to scalar masses in 
the ranges: 8 x 1(T 8 eV < m < 2.6 x 1(T 7 eV (for 1 + re = 10 2 ), 2.5 x 1(T 6 eV < m < 6 x 1(T 6 
eV (for 1 + re = 10 3 ), 8 x 1(T 5 eV < m < 2 x 1(T 4 eV (for 1 + re = 10 4 ), 2.5 x 1(T 3 eV 
< m < 6 x 1(T 3 eV (for 1 + re = 10 5 ), and 0.08 eV < m < 0.2 eV (for 1 + re = 10 6 ). Similarly, 
in figure 2, the resulting scalar mass range is: 1.4 x 10 _12 (1 + re) eV < m < 7 x 10~ 12 (1 + re) 
eV. A glance at these mass ranges 3 shows that the solutions in figure 1 and figure 2 exclude 
3 Scalar masses were obtained making [i e = 2 and p c = I0 8 g/cm 3 . For other choices, one should multiply 
the given numerical mass values by the factor (4p) 3 ( 10 a,gJ cm ;i ) 6 • 
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Figure 2: Equilibrium configurations as a function of the parameter £ L . As one moves along 
the curve from right to left, £l grows from 3 to £i. 

each other, i.e. either we have stars of type 1 or we have stars of type 2, but not both. 

Next issue on our agenda is to check whether the values for re and m make physical sense. 
As to re, we use the specific potential in equation (12) with a = rwyS" 2 ~ A-kGo 2 p c mR~ 2 
and use T c = 0(1) eV. We set the scale for the temperature T c at 1 eV because with this 
choice the temperature of the scalars will be below the cosmic background temperature at 
matter domination when formation of structure begins. This is reasonable since the scalars, 
which were in thermal equilibrium with photons in the early stages of the history of the 
Universe when temperatures were on the order of the hypothesized high energy unification 
scale, have not experienced the reheating associated to the various extinction processes (e.g., 
pp, e + e~) that have otherwise risen the temperature of photons. Little algebra and a 

short numerical calculation will convince us that g 2 varies in the interval 

''-^"--^'^fcw) 1 (36) 

for the values of re considered. Thus, an interaction which is extremely weak at a microscopic 
level is enhanced to considerable strength by the macroscopic coherence of the scalar medium. 

At this point one may wonder if it is realistic to think that the scalar medium can be kept 
at a temperature below leV during the formation of the white star while the temperature of 
the star itself can reach temperatures up to 10 4 — 10 5 eV. But, in fact, the scalars and the 
stellar matter are not in thermal contact. They do not share a common temperature. Indeed, 
the tiny coupling of scalars to nucleons implies an elastic scattering cross section off non 
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relativistic nucleons which is at most a^N ~ (9(lCT 64 cm 2 ). Clearly, a£ N is so much smaller 
than the nucleon spacing in a nucleus, that the nucleus is transparent and the cross section 
a^A on a nucleus with mass number A is ~ A a^. As a consequence, the largest possible 
collision rate of a scalar with nuclei 4 in stellar matter (avn) is so extremely small that not 
even in a hundred white dwarf lifetimes a scalar would have collided once with a nucleus. 
Needless to say, fixing the condensation temperature at about leV is by no means crucial. 
There is a wide band of temperatures for which the resulting scenario is phenomenologically 
sound. 

Also, the effects of this scalar force would go unnoticed by local gravity experiments. 
Indeed, the ratio of the vacuum potential (equation (13); we suppose that in a terrestrial 
environment, there is no appreciable amount of scalars trapped) to the newtonian potential 
is 

3 

Vvac - 6.9 x 1(T 23 ( ( -^-Y ( —) 2 e~ 2mr (37) 



V Newt ' V 10 35 7 V10-W/ \mr 

for r ^> m _1 , and 



for r <C vrC 1 . We will see below that the domain within which m _1 is contained is roughly 
~ 10~ 4 cm — 10 2 cm. No experiment searching for new forces has or will exclude in a foreseeable 
future such a tiny force [T3]. 

A comment on m comes next. The condensed phase consists of scalar particles with 
macroscopic de Broglie wavelength X^b ~ The Compton wavelength of the particles, 

on the other hand, should be much smaller than their de Broglie wavelength for the zero- 
point momentum to be negligible as compared to mass. In the low £i case this means: 
m ^> \ d g = 1CT 11 — 1CT 9 eV for (1 + k) between 10 2 and 10 6 , respectively. The masses 
associated with figure 1 show no conflict with this requirement, where even in the most 
unfavorable case (for (1 + k) — 10 2 ) m is almost a factor 10 4 bigger than the corresponding 
inverse de Broglie wavelength. In the high case one should have: m 3> X^ B = 10~ 12 — 10~ 10 
eV for (1 + k) between 10 2 and 10 6 , respectively. Obviously, the range 1.4 x 10 _12 (1 + k) 
eV < m < 7 x 10 _12 (1 + k) eV in figure 2, comfortably satisfies this requirement only 
for sufficiently high (1 + k). In short, probably the mass m can be reasonably constrained 
to lie between (9(1CT 7 — 1CT 6 ) eV and 0(0.1) eV. Or, in terms of Compton wavelength: 
~ 10~ 4 cm — 10 2 cm. Hence, the range of the 2-scalar force for masses such as above — if it 
were not for the condensation phenomenon — would be at most on the order of one meter. It 



4 With A ~ 20. 
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Figure 3: White dwarf masses and radii. Data points correspond to the white dwarfs Sirius 
B, Stein 2051, 40 Eri B and Procyon B in binaries, and to the white dwarfs in the Common 
Proper-Motion Pairs CD-38 10980, W485A, L268-92, G181-B5B and G156-64 [H]. Theoret- 
ical calculations lie within the two dotted lines. 

is only because the dispersion forces become of infinite range that they can coherently extend 
over a large piece of stellar matter. 

The values of m entertained in the preceding paragraphs are below or much below the 
temperature T ~ 1 eV. Hence the scalars are relativistic or ultrarelativistic as they should. 
Furthermore, T/m C Q as can be seen from: 

«^, 5xl o« (1 + K) l(^)^(^) ? (^) 4 (^) 3 W 

which follows from equation (5). 

We conclude with a summary. Two scalar exchange among nucleons produces spin inde- 
pendent Van der Waals type forces. In vacuum their range is about the Compton wavelength of 
the scalar particles exchanged. But when ordinary matter is embedded in a Bose-Einstein con- 
densate composed of the very same mediating scalars, those forces become infinitely ranged. 
The force acts then coherently over macroscopic extensions of bulk matter. The phenomenon 
arises as a combination of kinematics — 3-momentum exchange of the matter system with 
the scalar medium — and the collective effect of condensation. The main purpose of the 
present paper has been to find and study a physical system where this phenomenon of infinite 
widening of the range could be relevant. 

Scalars are fundamental ingredients of the SM of particle physics and completions thereof. 
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Figure 4: The two dotted upper lines limit the area where type 1 stars (for the choice of 
parameters: 1 + k — 10 2 and m = 1.5 x 10 -7 ) should lie. The two dotted lower lines 
correspond to the lines for actual white dwarfs in figure 3. 

Quite generally, the interactions of new matter with ordinary matter at the postulated new 
high energy unification scale in those beyond the SM scenarios will render low energy residual 
effects, however tiny, involving ordinary matter at rest. In particular, the dispersion forces 
due to the double exchange of scalars just mentioned will show up. These scalars, on the other 
hand, may constitute a component of dark matter. If this is the case, it might well be that 
chunks of scalar galactic dark matter in the condensed phase be trapped in stellar media. As 
a result, long range scalar forces would affect the equilibrium of those stars. Restricting our 
analysis to white dwarfs, we have explored the changes induced on their radii and maximal 
mass — i.e. the Chandrasekhar limit — by the presence of dispersion forces. Although figures 
1 and 2 show only this limit, in order to relate our results to actual white dwarf data, we do the 
following exercise. Figure 3 displays masses and radii for a handful of white dwarfs obtained 
from visual binaries or common proper-motion systems Also shown are two dotted 

lines that embrace the different model calculations (including the pioneering Chandrasekhar 
models) [13 EH . We can get a rough idea for where in the mass-radius plane our hypothetical 



M 



and 



on 



stars might lie, by scaling the band between dotted lines on figure 3 with l , ,,,,,, 
figures 1 and/or 2. As an example let us take 1 + k — 10 2 and m — 1.5 x 10 -7 eV which 



corresponds to a type 1 star with 



M 
M Ch 



0.815 and = 11 (see figure 1). Figure 4 shows 



the resulting band as well as, for comparison, the band for ordinary white dwarfs of figure 3. 
We see that, depending on the mass m of these scalars and strength nG of the potential, a 
distinctive feature of these star configurations would be that they should populate, relative to 
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ordinary white dwarfs, patently different regions on the (M, i?)-plane. A novel clue on dark 
matter might thus come from the identification of such anomalous white dwarfs. 
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